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Abstract. Theory is presented for the phase stability of mixtures containing 
nanospheres and non- adsorbing reversible supramolecular polymers. This was made 
possible by incorporating the depletion thickness and osmotic pressure of reversible 
supramolecular polymer chains into generalized free- volume theory, recently developed 
for investigating the phase behaviour of colloidal spheres mixed with interacting 
polymers [Adv. Colloid Interface Sci. 143 (2008) 1-47]. It follows that the fluidfluid 
phase stability region where reversible supramolecular polymer chains can be mixed 
with nanospheres is sensitive to the energy of scission between the monomers and to 
the nanoparticle radius. One can then expect the fluidfluid coexistence curves to have 
a strong dependence on temperature and that shifting of phase boundaries within 
a single experimental system should be possible by varying the temperature. The 
calculations reveal the width of the stability region to be rather small. This implies 
that phase homogeneity of product formulations containing reversible supramolecular 
polymers is only possible at low nanoparticle concentrations. 



PACS numbers: 81.16.Fg, 64.75.Xc, 82.70.Dd 



Phase stability of reversible polymers and nanospheres 2 
1. Introduction 

Reversible supramolecular polymers are polymer cliains tliat form due to spontaneous 
self-assembly of monomers [1]. Recent progress in supramolecular chemistry led 
to the development of several types of reversible supramolecular polymers including 
hydrogen-bonded equilibrium polymers [21 [3l HI El E] or coordination polymers that 
self-assemble via for instance metal ion-ligand group interactions [7J. Another class 
of such reversible polymers that behave similarly are worm-like surfactant micelles [8]. 
Such supramolecular polymers can be used to form a large diversity of self-organized 
multifunctional materials [9]. 




Figure 1. Sketch of the formation of reversible supramolecular polymers composed of 
individual monomers. By changing for instance temperature or monomer concentration 
the equilibrium constant attains a different value leading to a modified chain length 
distribution. 



The chain length (distribution) of these polymers is determined by the end-cap or 
scission energy and the polymer concentration [8]. The scission energy, U, is the energy 
required to create two new chain ends. A schematic picture of the formation of these 
reversible polymers is drawn in Fig. [TJ A larger monomer concentration and stronger 
attraction between the monomers enhance self-assembly into longer chains. Often the 
bonds form and break on time scales that are accessible experimentally. The fact 
that the bonds between the monomers are noncovalent provides solutions of reversible 
supramolecular polymers with a chain length distribution that depends on monomer 
concentration and temperature. Here a simple mean-field model [8] is adopted to 
describe the equilibrium chain length distribution. For sake of simplicity only athermal 
chains are considered. For a given segment volume fraction the segment volume 
fraction that has a chain length M, (pM, is then given by the exponential distribution: 

,„.^e-'/<«) . (1) 

The brackets (■) denote the ensemble average; (M) is the number- averaged chain length 
given by: 

(M) = ^ e^/^ , (2) 
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where U is the scission energy between the monomers (normahzed by kT) . Following van 
der Gucht et al. [TOj, the relative concentration of reversible supramolecular polymers 
(further denoted as equilibrium polymers) 0p is defined as 



with 0p = 1 defining the overlap concentration of equilibrium polymers at which the 
equilibrium polymer coils exactly fill the total volume. 

In this paper the effect of adding nanospheres to solutions of non-adsorbing 
equilibrium polymers is considered. Such nanoparticles can add a certain functionality 
to the supramolecular polymer solution. With the increasing interest in supramolecular 
polymers and widening possibilities to synthesize them, the phase stability of such 
mixtures will become increasingly relevant. As far as I am aware, only a single 
(preliminary) experimental study on the phase stability of nanoparticles mixed with 
(non-adsorbing) equilibrium polymers has been reported [llj. Knoben et al. examined 
2,4-bis(2-ethylhexylureido)toluene solutions mixed with stearylated silica nanospheres 
in cyclohexane and detected a liquid-liquid phase separation at equilibrium polymer 
concentrations at and above 10 g/L {ip > 0.01). 

The system of interest is sketched in Fig. |2l The equilibrium polymers are depleted 
from zones around the spheres termed depletion layers (dashed). In this paper, the 
mixture of nanospheres and equilibrium polymers in a solvent is studied using free 
volume theory. 



Figure 2. Schematic picture of nanospheres in a solution containing non-adsorbing 
equilibrium polymers solution. Note the broad chain length distribution of the 
equilibrium polymers. The effective depletion layers are indicated by the dashed circles. 
Upon overlap of depletion layers (the lower two spheres) the osmotic pressure acting 
upon the spheres is unbalanced. The resulting attractive force is indicated by arrows. 



(pp = ^ (M) 



(3) 
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In the early 1990s Lekkerkerker and others put forward free volume theory 
(FVT) for mixtures of spherical particles and (ideal) polymers |12l [13]. For solutions 
containing non-adsorbing polymers, the nanoparticles are surrounded by depletion layers 
[m [151 [IS]- Overlap of depletion layers results in attraction between the particles [17] 
and these attractive forces in turn induce phase transitions in nanoparticle plus polymer 
mixtures [HI [TH [201 Ell [H] ■ The phase diagram topology depends on the polymer-to- 
sphere size ratio. 

The power of FVT is that it is simple, insightful, and accurate (see for instance 
[22l [23l [211 [25]). Over the years, FVT has been shown to be flexible in that it can be 
used to describe the phase behaviour of mixtures of colloidal rods and polymers [26] , 
dispersions of colloidal spheres and (small) colloidal rods [27], charged spheres mixed 
with polymer chains [281 EH] , and mixtures of interacting polymers and colloidal spheres 
[301 [an [32]. 

Here I present generalized free volume theory for a mixture of spherical 
nanoparticles in a supramolecular polymer solution. The theory is outlined in section 
[21 and the results for the phase stability are presented in section [31 

2. Theory 



2.1. Semi-grand canonical potential 

Within generalized free- volume theory (GFVT), a nanoparticle-polymer mixture is 
described by the semi-grand-canonical ensemble. The system of interest is held in 
osmotic equilibrium with a reservoir (r). This reservoir is separated from the system 
via a hypothetical membrane that allows permeation of equilibrium polymers and 
solvent molecules but which does not allow the nanospheres to enter the reservoir. The 
equilibrium polymer concentration in the system is set by the chemical potential of the 
polymer chains in the reservoir, /ip. The system, with volume V at temperature T, 
consists of A^n nanospheres with radius a in a solvent with equilibrium polymers. The 
solvent can safely be considered as 'background' [121 [331 121] • An expression is needed for 
the semi-grand potential Q for the mixture of nanoparticles and equilibrium polymers in 
solution in order to be able to compute the thermodynamic properties of the mixture. At 
given system volume and temperature, the normalized semi-grand canonical potential, 
u = (3Qvn/V, can be written in normalized form as [301 EI] 

co = f + Up. (4) 

with 

where the reduced canonical free energy = PF^v^/V. Here 1//3 = k^T is the thermal 
energy, = Aira^/S is the volume of a nanosphere, and is the Helmholtz energy 
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of the pure nanoparticle dispersion. The nanoparticle concentration is expressed in the 
volume fraction rj = N^v^/V . 

In order to calculate the polymer concentration in the system one needs the free- 
volume fraction a for the polymer chains in the system, defined as 0(77, /xp = (Vfree)/^, 
where Vfree is the free-volume in the system not occupied by the nanospheres and their 
depletion zones. 

For the free volume fraction a the standard scaled-particle result is used [3U [35| 
ESI [131 SO]: 

a = (1 - r/) exp(-Ay - By^ - Cy^), (6) 

where y = r//(l - r]), A = {1 + 5/af - I, B = 3{6/af{6/a + 3/2), and C = 3{6/af. 
Here 6 is the thickness of the depletion zone around a nanosphere with radius a. From 
([5]) and ([6]) it follows the quantities S and dU/dcpp are needed to calculate cUp. These 
will be given below. 



2.2. Reversible supramolecular polymer contribution 

The ensemble-averaged chain length (M) as a function of the normalized polymer 
concentration 0p follows from inserting Q into (|2]): 

The size of the equilibrium polymers Re is defined as: 



It is noted that all length scales are normalized with the size of a polymer segment in 
this paper. 

The depletion thickness for dilute and semi-dilute equilibrium polymer solutions 
next to a fiat wall, 6p, was obtained previously [TO] , 

R, " 



'5p 



(9) 



In the dilute limit (0p — j- 0), 6p vanishes and in the semi-dilute limit (0p ^ 1), the 
depletion thickness follows the mean-field scaling relation 5p ~ 0p [TU] . 

It is important to also account for curvature effects to accurately describe the 
depletion thickness, 6, around a spherical particle. In the mean-field case 6 can be 
approximated accurately by the simple power law [21] : 

- = 0.842 f^V'' . (10) 
a \(^ / 

In Fig. [3] the characteristic polymer size [Re) and depletion thickness {5) of equilibrium 
polymers are plotted as a function of the relative polymer concentration (pp for a few 
values of the scission energy. It illustrates how the size of the equilibrium polymer 
chains (i?e) increases with 0p. The depletion thickness also increases (as it follows R^ 
in the dilute case) up to concentrations close to overlap (of order unity) above which 
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6 decreases; in the semi-dilute polymer concentration regime S follows the correlation 
length in a polymer solution |16]. As the scission energy is larger, the chains will be 
larger at given relative polymer concentration 0p. The same holds for the depletion 
thickness. 
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Figure 3. Concentration dependence of the depletion thickness (solid curve) and chain 
size (dotted curve) of the equilibrium polymers around a sphere for three values of the 
scission U energy as indicated. 



For the osmotic pressure 11 of equilibrium polymers the following mean-field result 
is used: [10] 

P = fiv,U = <P^ + ^^l . (11) 

Here P is the normalized osmotic pressure. The polymer concentration dependence of 
osmotic pressure is plotted in Fig. H] for various values of the scission energy U. The 
derivative of osmotic pressure with respect to polymer concentration is 

dP ^ dn f aV ^ ^ 

P^-T^ = hr [i + <^p] • (12) 



2.3. Free energy of a nanoparticle dispersion of hard spheres 

The canonical free energy of a fluid dispersion containing only hard spheres can be 
accurately described using: 

f = V HvA'M - V + ^^^^ , (13) 

where A = h/ y/ 2-71 m^kBT is the thermal wavelength, mn is the nanoparticle mass and h 
is Planck's constant. The first two terms on the right-hand side of this equation comprise 
the ideal contribution, whereas the last term describes the hard-sphere interaction, which 
follows the equation of state given by Carnahan and Starling [37] . 
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Figure 4. Osmotic pressure of equilibrium polymers for three values of the scission U 
energy as indicated. 



Standard thermodynamics enables computation of the osmotic pressure of the 
nanoparticle dispersion and the chemical potential of the nanoparticles. The chemical 
potential Jl^ = fi^/kT = df^/dr] follows as: 



=1^- + ^^^+ (1%)3 • 



(14) 



The dimensionless pressure P° = P^VrJkT can subsequently be computed from P° 
- /°, giving 



P 

the well-known Carnahan-Starling pressure. 



(15) 



2.4. Fluid- fluid binodals 

All ingredients are now available to calculate the fluid-fluid phase behaviour of the 
nanoparticle-equilibrium polymer mixture by solving the coexistence equations 

V'i = J^n , (16) 

and 



pi ^pii 

tot tot 



(17) 

for the chemical potential of the nanospheres and the total osmotic pressure of a phase 
I in equilibrium with a phase /J. 

The phase diagrams for hard spheres plus interacting polymers using the general 
expression (jl]) and its ingredients now follow from computing the chemical potential 
/In = {du/dri) and total pressure Ptot = vl'-n — ^, 



^ .0 / " I 9P' 

/^n = + / 9 

Iq 



R' 

R' I ^'^P 



Phase stability of reversible polymers and nanospheres 



Ptot = P'+ I h 



kR' 



8 

(19) 



with g and h given by: 

g = e-<^ {l + [l + y]{A + 2By + 3Cy^)} (20) 

and 

h = e-^ {l + Ay + 2By'^ + 3Cy^} . (21) 

with A, B and C defined in (and below) equation (|6]) and where Q is given by 

Q = Ay + By^ + Cy\ (22) 

The normahzed osmotic pressure of the equihbrium polymers in the reservoir is given by 
Eq. f|TT]) and the osmotic compressibility by Eq. f|T2|) . Coexistence curves then follow 
from Eq. f|T6l) and Eq. f|T7j) by using the common tangent construction. 
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Figure 5. Diniensionless semi-grand potential a; as a function of nanoparticle volume 
fraction rj. The common tangent construction identifies the binodal compositions rj^ 
and 77^^. In this way one can determine the phase coexistence in mixtures of hard 
nanospheres and equilibrium polymers. 

In Fig. |5]it is demonstrated how the (normalized) semi-grand potential depends on 
nanosphere particle volume fraction under conditions where the fiuid phase exhibits an 
instability as follows from the local maximum in u. Such a curve depends on the reservoir 
equilibrium polymer concentration, the nanoparticle size and on the scission energy. The 
common tangent construction in Fig. |5] allows to determine the compositions of the two 
coexisting phases. From these binodals the full phase diagram can be constructed. 
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Figure 6. Fluid- fluid coexistence curves of equilibrium polymers mixed with 
nanospheres at fixed scission energies U = 20 for various values of the nanoparticle 
radius a —5 (sohd curve), a = 10 (dashed curve) and a = 20 (dotted curve). Dots are 
the critical points. 

3. Results 

In Fig. |6] fluid- fluid binodal curves are presented for nanosplieres of a = 5, 10 and 20 
witfi sizes given in units of tlie equilibrium polymer monomer size. The scission energy 
here was fixed at 20 (kT). The left panel represents the relative polymer concentration 
in the reservoir (0p on the ordinate versus the nanoparticle volume fraction (77) on the 
abscissa. In this representation, it follows that relatively larger polymer chains (smaller 
a) shift the fluid-fluid binodal curves to higher relative polymer concentrations. It is 
clear that the binodal sensitively depends on a. The relative polymer concentration in 
the system, 0p, follows from (pp = acpp. The resulting binodal curves are quite close to 
each other (right panel in Fig. [6]). This is due to the fact that the free volume fraction 
depends on the relative polymer concentration and passes through a maximum around 
(j)p = 0.6. This result shows that the stability region is rather limited for mixtures of 
nanospheres and equilibrium polymers. 

In Fig. [7] fluid- fluid binodal curves are shown for mixtures of nanospheres with fixed 
size a = 10 mixed with equilibrium polymers for various values of the scission energy 
U, the change of which is experimentally accessible via manipulation of temperature. 
It follows that increasing U leads to a significant shift upwards of the fluid-fluid 
binodals when the relative equilibrium polymer concentration is plotted in the reservoir 
representation (left panel of Fig. [7]). A larger value for U leads to a larger average 
equilibrium polymer size. Therefore these results are in agreement with (G)FVT results 
for (normal) polymer-sphere mixtures |13l UHl |30] . Conversion towards system polymer 
concentrations as in practice (right panel of Fig. [7]) shows the stability region is small 
and the binodal curves are close. 

What is observed in practice? The relative polymer concentration is insightful since 
it identifies the concentration regime; dilute regime for (pp < 1, semi-dilute regime if 0p 
> 1. The quantity (pp is very sensitive to the scission energy U, while the segment 
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Figure 7. Fluid-fluid coexistence curves of equilibrium polymers mixed with 
nanospheres of radius a = 10 for scission energies U = 15 (dotted curve), U = 20 
(dashed curve) and U — 25 (solid curve). Dots: critical points. 



concentration ip is of course unaffected. The phase diagram of Fig. [7] (right panel) is 
now plotted in terms of the segment concentration on the ordinate, as represented in 
Fig. [HI It follows that variations of U can have a dramatic effect on the binodal fiuid- 
fiuid coexistence curves observed in experiment. Changing the relative size of polymer 
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Figure 8. Fluid-fluid coexistence curves of equilibrium polymers in system segment 
concentrations ip mixed with nanospheres of radius a for scission energies U — 15 
(dotted curve), U — 20 (dashed curve) and U — 25 (solid curve). Dots are the critical 
points. 

and spherical particle for mixtures of (normal) polymers mixed with nanoparticles 
requires significant experimental effort; one needs a different system with either particles 
or polymers of another size. This differs from changing U in mixtures of reversible 
supramolecular polymers and nanospheres in the sense that one only needs to modify 
the temperature. With a single mixture one can then scan through the phase diagram 
by diluting, concentrating and changing the temperature. No quantitative comparison 
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can yet be performed with experimental work. The single experimental study that 
is available pT] used chain stoppers, which means we can not directly compare with 
their system. Further, their observations of (in)stability were quite limited. It will be 
interesting to compare the theory presented in this work to future phase diagrams on 
mixtures of nanoparticles and equilibrium polymer chains. 

No remark was made yet on the critical point. For short-range attractions it is 
known that mean-field theories incorrectly predict the critical region, whereas long- 
ranged attractions are well described. In Figures [6] and [7] it can be observed that 
increasing U and decreasing a imposes similar effects on the critical points and the 
binodals. These critical (pp and rj values are plotted as a function of a and f/~^ in a 
single plot in Fig. [91 It follows that, when scaled properly, the influence of a and 
are very similar and that one can modify the binodals in a similar fashion via adjusting 
a or U. 

a 

5 10 15 20 25 
0.5-r-, ■ 1 ■ 1 ■ 1 ■ 1 0.5 




Figure 9. Critical curves for 0p and as a function of a and 1/U demonstrating their 
effects lead to similar shifts of the fluid-fluid binodal curves. 



4. Concluding remarks 

Theory was presented for the fluid-fluid phase transition of mixtures of dispersed 
nanospheres in solution containing equilibrium polymers. This was based upon 
generalized free volume theory originally developed for spherical particles mixed with 
interacting monodisperse polymer chains. The equilibrium polymers in solution were 
treated as athermal chains using a mean-field model. Within this model analytical 
expressions are available for the size, chain length (distribution), osmotic pressure 
and depletion thickness. In the dilute monomer concentration regime, the depletion 
thickness increases with monomer concentration while depletion thickness decreases in 
the semi-dilute equilibrium polymer concentration regime. A key parameter that sets 
the averaged chain length of the polymers is the scission energy; the energy required to 
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create two new chain ends. In practice the scission energy can be varied by adjusting the 
temperature. An increase in scission energy results in increases in depiction thickness, 
while the osmotic pressure decreases at given relative polymer concentration. 

The generated fluid-fluid binodal curves can partly be interpreted using previous 
findings from (generalized) free volume theory. These results can provide new insight 
into experimental observations. Upon varying the scission energy via the temperature 
one can, in combination with dilution and concentration of the dispersion, scan through 
the scission energy-dependent phase diagram. Finally, it was found that the critical 
point scales in a similar way with scission energy as with the inverse of the nanosphere 
radius. 

In a later contribution, the effect of polymer solvency, as well as an account of the 
fluid-solid equilibrium, is planned. 
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